DESY 98-025 
March 1998 



o 



^1" 



Hadron Spectra and Spectral Moments in the Decay 

B X s £ + £~ using HQET 



A. Ali* and G. HiW 

00 ■ Deutsches Elektronen-Synchrotron DESY, Hamburg 

ON 

S3 



Abstract 

We compute the leading order (in a s ) perturbative QCD and power (l/m|) corrections to the 
hadronic invariant mass and hadron energy spectra in the decay B —> X s £ + £~ in standard model 
using the heavy quark expansion technique (HQET). Results for the hrst two hadronic moments 
' (^h) ano - C^h)' n — 1' 2, are presented here working out their sensitivity on the HQET parameters 

Ai and A. Data from the forthcoming B facilities can be used to measure the short-distance 



contribution in B — > X s £ + £ and determine the HQET parameters from the moments (S^). This 
can be combined with the analysis of semileptonic decays B — > X£vg to determine them precisely. 
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The semileptonic inclusive decays B — ► X s £ + £~ , where I = e , fj, ,t , offer, together with the 
radiative electromagnetic penguin decay B — > X s + 7, presently the most popular testing grounds for 
the standard model (SM) in the flavor sector. Concentrating on the decay B — > X s £ + £~, we recall 
that the first theoretical calculations including partial leading order QCD corrections were reported 
a decade ago in emphasizing the sensitivity of the dilepton mass spectrum and decay rates to 

the top quark mass in the short-distance contribution. Since the pioneering papers [|]-[|, a lot of 
theoretical work has been done on the decay B — ► X s £ + £~. In particular, the complete leading order 
perturbative corrections in the QCD coupling constant a s to the dilepton invariant mass spectrum 
and so-called forward-backward (FB) asymmetry of the leptons [||, have been calculated in refs. ||||] 
and 0], respectively. In addition, leading order power corrections in to the decay rate, dilepton 

invariant mass spectrum and the FB asymmetry have been studied Q, using the heavy quark expansion 
technique (HQET), correcting an earlier derivation of the dilepton spectrum in this decay H|. The 
power corrected dilepton mass spectrum and FB asymmetry have also been derived for the massless 
s-quark case recently in Q, confirming the results in |7|]. Corrections of order l/m 2 to the dilepton 
mass spectrum away from the resonant regions have also been worked out [p~0|,p~T|] . 

In this letter, we study the hadron spectra and moments in B — > using HQET. The study of 



the decay B — > Xivi in this context has received a lot of interest [12-17] . The hadronic invariant mass 
spectra in B — ► X s £ + £~and B — > X u £vi have striking similarities and differences. For example, both of 
these processes have at the parton level a delta function behavior dT/dso oc S(sq — rriq), q = u,s, where 
so is the partonic invariant mass. Thus, the entire invariant mass spectrum away from so = mi is 
generated perturbatively (by gluon bremsstrahlung) and through S-hadron non-perturbative effects. 
The latter are already present in 0(a®), as can be seen by the relation between the b quark mass and 
the B meson mass. In HQET this takes the form = nib + A — (Ai + 3A2) /2mi, + .... The quantities 



Ai, A2 and A are the HQET parameters [18, 19, a]. Keeping, for the sake of simplicity just the A term, 
the hadronic invariant mass Sh is related to sq and the partonic energy E by S H = s + 2AE + A 2 . 
This gives rise to a non-trivial spectrum in the entire region A 2 < Sh < M B . Hence, measurements 
of these spectra would lead to direct information on the QCD dynamics and a better determination 
of the non-perturbative parameters, such as the HQET parameters A and Ai- Following this line of 
argument, the sensitivity of the lepton and hadron energy spectra on these parameters in the decays 



B — > X£i/£ has been studied quantitatively in literature [12,13,16] ; Photon energy moments in the 
decay B — > X s + 7 have also been worked out in p0[| . Present status of the HQET parameters is 
reviewed in [21]. There is a fair amount of theoretical dispersion on A and Ai and it will be very 
instructive to get independent and complementary information on these parameters from other B 
decays. 

We report here a calculation of the hadron spectra in the decay B — > X s £ + l~ . Leading order 
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(in a s ) perturbative QCD and power (1/m?) corrections to the hadronic invariant mass and hadron 



"bt 
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energy spectra in this decay are computed at the parton level. Including both the 0(l/m£) and 
0(a s ) terms generates hadron spectra with contributions of 0(A/ms), C(a s A/ms), 0(Xi/m B ) and 



0(\2/m 2 B ). Relegating the detailed hadronic profile to a subsequent publication [22], here the power- 
and perturbatively corrected hadronic spectral moments {S H ) and {E H ) are presented for the first two 
moments n = 1, 2. The former are sensitive to the HQET parameters A and Ai and we work out this 
dependence numerically, showing that these moments would provide an independent determination of 
the HQET parameters in B —* X s £ + £~ . The theoretically constrained contours in the (A — Ai) plane 
in B — > X s £ + £~&re compared with the corresponding one from an analysis of the power corrected 
lepton energy spectrum in B —* Xlv^ fl3]] . We argue that a simultaneous analysis of the moments 
and spectra in B —* X s £ + £~ ?m<\ B — > X£v% will allow to determine the HQET parameters with a high 
precision. 

We start with the definition of the kinematics of the decay at the parton level, b(pb) — * s (Ps)(+g(Pg))+ 
+ 1~ (p-), where g denotes a gluon from the 0(a s ) correction. The corresponding kinematics 
at the hadron level can be written as: B(ps) — ► X s (j>h) + £ + {p+) + t~{pJ). We define by q and s 
the momentum transfer to the lepton pair and the invariant mass of the dilepton system, respectively, 
q = p + +P-, s = q 2 ; the dimensionless variables with a hat are related to the variables with dimension 
by the scale m&, the 6-quark mass, e.g., s = rh s = ^ etc.. Further, we define a 4-vector v, 
which denotes the velocity of both the 6-quark and the i?-meson, pb = rribV and ps = tubv. We 
shall also need a variable u, which is defined as u = —{pb — P+) 2 + (Pb ~ P-) 2 ; with the scaled vari- 
able u = u/m 2 satisfying the kinematic relation u = 2v • (p+ — p~). The hadronic invariant mass 
is denoted by Sh = p\ and Eh denotes the hadron energy in the final state. The corresponding 
quantities at parton level are the invariant mass sq and the scaled parton energy xq = without 
gluon bremsstrahlung this simplifies to sq = m 2 and xq becomes directly related to the dilepton in- 
variant mass, xq = 1/2(1 — s + rh 2 ). From momentum conservation, the following equalities hold in 
the 6-quark, equivalently £>-meson, rest frame (v = (1,0,0,0)): 

xo = 1 — v ■ q , so = 1 — 2v ■ q + s , 
Eh = rriB — v-q, Sh = m 2 B - 2m B v ■ q + s . (1) 

The relation between the kinematic variables of the parton model and the hadronic states is, using 
the HQET mass relation mg = nib + A — l/2m;,(Ai + 3A2) + • ■ ., given as 

p t A1 + 3A2 / -. A! + 3A 2 ^ 

Eh = A h tub — Ah \ x + ... , 

IrriB \ 2m b J 

S H = m 2 + A 2 + (m 2 B - 2Am B + A 2 + Ai + 3A 2 ) (s - m 2 ) 

+ (2Am B - 2A 2 - Ai - 3A 2 )a;o + . . . , (2) 



where the ellipses denote terms higher order in l/m&. The quantity A2 is known precisely from the 
B* — B mass difference, with A2 — 0.12 GeV 2 . The other two parameters are of interest here. 

The effective Hamiltonian governing the decay B — > X s £ + £~ , obtained by integrating out the top 
quark and the bosons, is given as 

v ts hb 

,i=l 
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(3) 



where L and R denote chiral projections, L(R) = 1/2(1 + 75), Vij are the CKM matrix elements and 
the CKM unitarity has been used in factoring out the product V t * s Vtb- The operator basis is taken 
from where also the Four-Fermi operators 0±, . . . , 0% and the chromo-magnetic operator 0% can 
be seen. Note that 0% does not contribute to the decay B — > X s £ + £~in the approximation which 
we use here. The Ci(n) are the Wilson coefficients, which depend, in general, on the renormalization 
scale except for C\q. Their numerical values are given in Table |. 



Ci 


c 2 


C 3 


<?4 


c 5 


c 6 




c 9 


C10 


-0.240 


+1.103 


+0.011 


-0.025 


+0.007 


-0.030 


-0.311 


+4.153 


-4.546 



Table 1: Values of the Wilson coefficients used in the numerical calculations corresponding to the 
central values of the parameters given in eq. Here, C^ 31 = C 7 - C 5 /3 - C 6 , and for C 9 we use 

the NDR scheme. 



With the help of the above expressions, one can express the Dalitz distribution in B — > X s 
dr 1 Gp 2 a 2 m b A lrr>frr l2 „ T ___ (mh tLW , mR t rv 



b —\v*v tb \ 2 2im(rV^ + rV^ 



as: 



(4) 



dudsd(v-q) 2 tub 2tt 2 256 

where the hadronic and leptonic tensors T L / R ^ V and L L I R ^ U are given in |7|]. Using Lorentz decom- 
position, the tensor T„„ can be expanded in terms of three structure functions Tj, 

-T\ g^ u + T 2 v v + T 3 v a q 13 



'- )1U 



(5) 



where the ones which do not contribute to the amplitude in the limit of massless leptons have been 
neglected. 

Concerning the 0(a s ) corrections to the hadron spectra, we note that only O9 = e 2 /(167r 2 )s Q 7 M L& Q £7^ 
is subject to such corrections. These can be obtained by using the existing results in the literature as 
follows: The vector current O9 can be decomposed as V = (V — A)/2 + (V + A)/2. Note that the 
(V — A) and (V + A) currents yield the same hadron energy spectrum |^] and there is no interference 
term present in this spectrum for massless leptons. So, the correction for the vector current case can 



be taken from the corresponding result for the charged (V — A) case [24 

3 



We have calculated the order a s perturbative QCD correction for the hadronic invariant mass in 
the range m 2 < so < 1. Since the decay b — ► s + £ + + i~ contributes in the parton model only at 
so = m 1i only the bremsstrahlung graphs b — ► s + g + ^ + + ^~ contribute in this range. This makes 
the calculation much simpler than in the full so range including virtual gluon diagrams. Also for this 
distribution, the results can be taken from the existing literature. As the starting point, we use the 
Sudakov exponentiated double differential decay rate g^j^, derived for the decay B — ► X u lvi in [fPfl , 
which we have checked, after changing the normalization for B — ► X s £ + £~. Defining the kinematic 
variables (x,y) as q 2 = x 2 m 2 , v ■ q = (x + 1/2(1 — x) 2 y)m;,, the Sudakov-improved Dalitz distribution 
is given by 

^-(B^X s £+r) = _^x(l-x 2 ) 2 (l + 2x 2 )exp(-^ln 2 (l- y )) (6) 

f 4a s ln(l - y) r 2a i 2a s dff \ 2 



where the functions G{x) and -ff(?/) can be seen in [14]. The constant Bq is defined below. 

The most significant effect of the bound state is the difference between mg and m& which is 
dominated by A. The spectrum is obtained after changing variables from (x,y) to (q 2 ,Sn) 
and performing an integration over q 2 . It is valid in the region wg 1 " 8 ^!" 18 < Sh < (or 
m^A < < neglecting m s ) which excludes the zeroth order and virtual gluon kinematics 

(so = m ^)- The hadronic invariant mass spectrum thus found depends rather sensitively on mj (or 
equivalently A). An analogous analysis for the decay B — > X u £i/£ has been performed in [15]. 

The hadronic tensor in eq. (||) can be expanded in inverse powers of nib with the help of the HQET 
techniques P Jl^Jl9t[ . The leading term in this expansion, i.e., 0(m®) reproduces the parton model 
result. In HQET, the next to leading power corrections are parameterized in terms of Ai and A 2 . The 
contributions of the power corrections to the structure functions Tj has been calculated up to (but not 
including) 0(l/mf) and given in 0. After contracting the hadronic and leptonic tensors and with the 
help of the kinematic identities given in eq. (jl]), we can make the dependence on x$ and so explicit, 

I^V L L / R " U = m b 2 { 2 (1 - 2xo + h)Ti L ' R + [x 2 - \u 2 - so] T 2 L ' R T (1 - 2x + h)u T 3 ^} 

and with this we are able to derive the double differential power corrected spectrum da ,^g for B — ► 



X s £ + £ . Integrating eq. (||) over u first, where the variable u is bounded by — 2wxq — sq < u < 



g — so, we arrive at the following expression [22]: 

— -^BolmJx 2 - s I (1 - 2x + Sq)T x (s , x ) + ^^T 2 (s , x ) 1 + C(Aja s ) . (7) 



dxo dso vr v I 3 

As the structure function T 3 does not contribute to the branching ratio, we did not consider it in 
our present work. The functions Ti(so,xo) and ?2(so,xo) have been derived by us after a lengthy 
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calculation and the resulting expressions are too long to be given here. They can be seen together 
with other details of the calculations in [J22[] . 

The branching ratio for B — * X s £ + £~is usually expressed in terms of the measured semileptonic 
branching ratio B s i for the decay B — > X c £u£. This fixes the normalization constant Bq to be, 

R _ R 3» 2 \vsy*\ 2 i (R) 

°°-^167r2 \ Vcb f /(m c ) K (m c )' W 
where f(m c ) is the phase space factor for T(B — ► X c £vi) and the function K(m c ) accounts for both 
the 0(a s ) QCD correction to the semileptonic decay width |2(| and the leading order (1/m;,) 2 power 
correction |jl8| 1. They are given explicitly in fl2|] . 

The hadron energy spectrum can now be obtained by integrating over sq. The kinematic boundaries 
are given as: max(rh 2 , —1 + 2xo + 4mf) < sq < Xq, rh s < xq < 3(1 + JTt 2 , — 4m 2 ). Here we keep m\ 
as a regulator wherever it is necessary. Including the leading power corrections, the hadron energy 



spectrum in the decay B — > X s ^ + ^~is derived by us and given in 22]. 

The lowest spectral moments in the decay B — » X s ^ + ^~at the parton level are worked out by 
taking into account the two types of corrections discussed earlier, namely the leading power l/m& and 
the perturbative 0(a s ) corrections. To that end, we define: 



ux- 

for integers n and m. These moments are related to the corresponding moments (xq l (sq — m 2 s ) n ) 
obtained at the parton level by a scaling factor which yields the corrected branching ratio B = 
BoM^'p. Thus, (xtf{§ — m 2 ) n ) — ^ We remind that one has to Taylor expand the 

correction factor Bq/B in terms of the 0(a s ) and power corrections. The moments can be expressed 
as double expansion in 0(a s ) and l/m& and to the accuracy of our calculations they can be represented 
in the following form: 

= D^ m) + ^C 9 2 A^ + Ai£>J n,m) + \ 2 D^ m) , (10) 
with a further decomposition into pieces from different Wilson coefficients for % = 0, 1, 2: 

D {n,m) = a in, m ) c eS 2 + ^ + ^) C f + (U) 

The terms <-y^ n ' m ) and fi^ 1 '" 1 ^ in eq. ( JXT] ) result from the terms proportional to Re{C^)C^ and |Cg^| 2 
in eq. ([/]), respectively. The explicit expressions for a [ n ' m \p^ m) 

The leading perturbative contributions for the hadronic invariant mass and hadron energy moments 
can be obtained analytically, 

A (o,o) = 25 ~ 4 ^ 2 4(1,0) = 4(2,0) = J_ 

9 675 ' 486 ' 

(0>1) 1381 - 210^ 2 (0 , 2) _ 2257 - 320^ 2 

A ~ 1350 ,A " 5400 ' [U) 
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The zeroth moment n = m = is needed for the normalization and we recall that the result for A^ ' ^ 
was derived by Cabibbo and Maiani in the context of the 0{a s ) correction to the semileptonic decay 
rate B — » Xivg quite some time ago f p6fl . Likewise, the first mixed moment A^ 1 ' 1 ' can be extracted 
from the results given in [12| for the decay B — > Xlv^ after changing the nor malization, AW = 3/50. 



For the lowest order parton model contribution D^ 1 ' 171 ^, we find, in agreement with [12], that the first 
two hadronic invariant mass moments (so — ml), ((so— ml) 2 ) and the first mixed moment (xo(so—m 2 )) 
vanish: D ( Q n ' 0) = 0, for n = 1, 2 and D { Q X,1) = . 

We can eliminate the hidden dependence on the non-perturbative parameters resulting from the 
6-quark mass in the moments A^+J- with the help of the HQET mass relation. As m s is of order 



AqcD) to be consistent we keep only terms up to order mj/m 2 . P7|| . An additional m^-dependence is 
in the mass ratios mi = With this we obtain the moments for the physical quantities valid up to 
0(a s /m B , l/m B ), where the second equation corresponds to a further use of m s = O(Aqcd)- We get 
for the first two hadronic invariant mass moments [] 

(S H ) = m 2 s + A 2 + (m| - 2Am B ) (s - m 2 s ) + (2Am B - 2A 2 - Ai - 3A 2 )(^o) , 

(S 2 H ) = m 4 s + 2A 2 m 2 + 2m 2 (m 2 B -2Am B )(so-m 2 ) + 2m 2 (2Am B -2A 2 -X 1 -3X 2 )(xo) 

+ [m% - 4Am B )((s - m 2 ) 2 ) + 4A 2 m|(x 2 ) + 4Am|(x (s - m 2 )) , (13) 
= (m% - 4Am|)((s - m 2 ) 2 ) + 4A 2 m|(^) + 4Am|(x (so - m 2 s )) , 

and for the hadron energy moments: 

(Eh) = A-±±^ + (m B -A + ±±^-)( X0 ), 
2m B V 2m B J 

(E 2 H ) = A 2 + (2Am B -2A 2 -Ai-3A 2 )(xo) (14) 
+(m| - 2Ami3 + A 2 + Ai + 3A 2 )(a?g) . 

One sees that there are linear power corrections, 0(A/m B ), present in all these hadronic quantities 
except (Sff) which starts in ^f^^- 



Using the expressions for the HQET moments derived by us [22], we present the numerical re- 
sults for the hadronic moments in B — > X s £ + l~ . The parameters used in arriving at the numerical 
coefficients are given in Table |l| and specified below: 

m w = 80.26 (GeV), m z = 91.19 (GeV), sin 2 9 W = 0.2325 , 
m s = 0.2 (GeV), m c = 1.4 (GeV), m 6 = 4.8(GeV), m t = 175 ± 5 (GeV), 
H = m b +™ b b/2 , a' 1 = 129, a s (m z ) = 0.117 ± 0.005 , B si = (10.4 ± 0.4)%. (15) 



x Our first expression for (Sfj), eq. does not agree in the coefficient of (so — m 2 s ) with the one given in (their 
eq. (4.1)). We point out that m% should have been replaced by m| in this expression. This has been confirmed by Adam 
Falk (private communication). Dropping the higher order terms given in their expressions, the hadronic moments in 
HQET derived here and in Q agree. 
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Inserting the expressions for the moments calculated at the partonic level into eqs. ( |i~3| ) and (p+|), we 

find the following expressions for the short-distance hadronic moments, valid up to 0(a s /m B , l/m%): 

, „ , <, ,m 2 e a, A a, A A 2 Ai Ao , 

(S H ) = m|(— f +0.093— -0.069 -+0.735 + 0.243^- + 0.273^- - 0.513—=-) , 

m B 7T TUB 7T Dig 1Tl B 771 B 771 B 

(Si) = mi (0.0071— +0.138— — + 0.587^- -0.196^-), (16) 
7r ms vr 

(E H ) = 0.367m B (l + 0.148— -0.352— — + 1.691— +0.012^ +0.024^i- + 1.070^§-), 

7r mg 7r mfl 

, o , 9 . a, A a, A A 2 Ai Aj , 

(E| = 0.147m|(l + 0.324— -0.128 - + 2.954 + 2.740^- - 0.299^- + 0.162—=-) . 

7r i m,B m B m B m B 

Concerning the non-perturbative parts related to the cc loop in B — > X s £ + ^~, it has been suggested in 
[ |l0| ] that an 0(AQ CD /m 2 ) expansion in the context of HQET can be carried out to take into account 
such effects in the invariant mass spectrum away from the resonances. Using the expressions (obtained 
with m s = 0) for the 1/m 2 amplitude, we have calculated the partonic energy moments A(xq), which 
correct the short-distance result at order \2jrr? c : 



B 256C9A9 /•V2(l-4mf) 



F W (cfff(3-2xo) + 2C?ff^±^±M 



where r = (1 — 2xo)/4m 2 and F(r) is given in [10]. The invariant mass and mixed moments give 



zero contribution in the order we are working for m s = 0. Thus, the correction to the hadronic mass 
moments are vanishing, if we further neglect terms proportional to — and -^-Aj, with i = 1,2. For 
the hadron energy moments we obtain numerically 



A(4)i« = m|A(x 2 ) = -0.013 GeV 2 , (17) 



A(E H ) 1/m 2 = m B A(x ) = -0.007 GeV , 

>l/mi 

leading to a correction of order —0.3% to the short-distance values presented in Table [2]. 

With the help of the expressions given above, we have calculated numerically the hadronic moments 
in HQET for the decay B — * X s £ + £~ , £ = /i, e and have estimated the errors by varying the parameters 
within their ±ler ranges given in eq. (|T5|). They are presented in Table [2] where we have used A = 
0.39 GeV and Ai = -0.2 GeV 2 . Further, using a s (m b ) = 0.21 and A 2 = 0.12 GeV 2 , the explicit 
dependencies of the hadronic moments given in eq. ( |T6| ) on the HQET parameters Ai and A can be 
worked out. 

A A 2 A 

(S H ) = 0.00557711(1 + 132.61 + 44.14^+49.66^-), 

tub m B m B 

A A 2 A 

(S H ) = 0.00048mfj(l + 19.41 + 1223.41^- - 408.39^-) . (18) 

ms m B m B 

As expected, the dependence of the energy moments (E H ) on A and Ai is very weak and we do 

not show these here. While interpreting these numbers, one should bear in mind that there are 
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two comparable expansion parameters A/ mg and a s /n and we have fixed the latter in showing the 
numbers. The correlations on the HQET parameters Ai and A which follow from (assumed) fixed 
values of the hadronic invariant mass moments (Sh) and (S H ) ar e shown in Fig. |]. We have taken 
the values for the decay B — > X s /j + \i~ from Table [2] for the sake of illustration and have also shown 
the presently irreducible theoretical errors on these moments following from the input parameters mt, 
a s and the scale fx, given in eq. (|l5|). The errors were calculated by varying these parameters in the 
indicated range, one at a time, and adding the individual errors in quadrature. As the entries in 
Table ^ are calculated for the best-fit values of Ai and A taken from the analysis of Gremm et al. 

for the electron energy spectrum in B — > Xiujg, there is no surprise that these curves meet at 
this point. This exercise has to be repeated with real data in B — ► X s £ + £~to draw any quantitative 
conclusions. Using the CLEO cuts on hadronic and dileptonic masses [pS||, we estimate that 0(200) 
B — ► X s £ + £~(£ = e,n) events will be available per 10 7 BB hadrons f22]| . So, there will be plenty of 
B — > X s £ + £~ decays in the forthcoming B facilities to measure the correlation shown in Fig. |l|. 

The theoretical stability of the moments has to be checked against higher order corrections and the 
error estimates presented here will have to be improved. The "BLM-enhanced" twodoop corrections 

proportional to a 2 /3o> where /3q = 11 — 2raj/3 is the first term in the QCD beta function, can be 



included at the parton level as has been done in other decays [12.3C] , but not being crucial to our point 
we have not done this. More importantly, higher order corrections in a s and are not included 

here. While we do not think that the higher orders in a s will have a significant influence, the second 
moment (Sjj) is susceptible to the presence of 1/mf corrections as shown for the decay B — > Xlv$ 
[p]]. This will considerably enlarge the theoretical error represented by the dashed band for (S#) in 
Fig. |l|. Fortunately, the coefficient of the A/ms term in (Sh) is large. Hence, a good measurement of 
this moment alone constrains A effectively. Of course, the utility of the hadronic moments calculated 
above is only in conjunction with the experimental cuts. Since the optimal experimental cuts in 
B — > X s £ + £~ remain to be defined, we hope to return to this and related issue of doing an improved 
theoretical error estimate in a future publication. The power corrections presented here in the hadron 
spectrum and hadronic spectral moments in B — > X s £ + £~&re the first results in this decay. 



HQET 


(Sh) 


(Sh) 


(Eh) 


(El) 




(GeV 2 ) 


(GeV 4 ) 


(GeV) 


(GeV 2 ) 


e + e~ 


1.64 ±0.06 
1.79 ±0.07 


4.48 ± 0.29 
4.98 ±0.29 


2.21 ±0.04 
2.41 ± 0.06 


5.14 ±0.16 
6.09 ±0.29 



Table 2: Hadronic spectral moments for B — > X s /j> + and B — » X s e + e~ in HQET with A = 0.39 GeV 
and Ai = — 0.2 GeV 2 . The quoted error results from varying [A,a s and the top mass within the ranges 
given in eq. (172 
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Figure 1: (Sh) (solid bands) and (S H ) (dashed bands) correlation in (X\-A) space for (Sh) = 
1.64 db 0.06 GeV 2 and (S H ) = 4.48 ±0.29 GeV 4 , corresponding to the values in Table |j. The curves 
are forced to meet at the point Ai = —0.2 GeV 2 and A = 0.39 GeV. The correlation from the analysis 
of the decay B — » X£vg from ref. [F3jJ is also shown here (ellipse). 



In summary, we have calculated the 0(a s ) perturbative QCD and leading 0(l/mb) corrections to 
the hadron spectra in the decay B — ► X s £ + £~ , including the Sudakov-improvements in the perturbative 
part. Hadronic invariant mass spectrum is calculable in HQET over a limited range Sh > m^A and 
it depends sensitively on the parameter A (equivalently m^). These features are qualitatively very 
similar to the ones found for the hadronic invariant mass spectrum in the decay B — * X u £vi [|D| . 
The l/m^-corrections to the parton model hadron energy spectrum in B — > X s £ + £~&ie small over 
most part of this spectrum. However, heavy quark expansion breaks down near the low end-point of 
this spectrum and near the cc threshold. We have calculated the spectral hadronic moments (E' H ) 
and (S H ) for n = 1,2 and have worked out their dependence on the HQET parameters A and Ax- 
The correlations in B —* X s £ + l~are shown to be different than the ones in the semileptonic decay 
B — > X£u£. This allows, in principle, a method to determine them from data in B — » X s £ + £~. We 
show the kind of constraints following from a Gedanken experiment in B — > X s £ + l~&nd the present 
analysis of data in B — > X£f£ [13] to illustrate this point. 

We thank Christoph Greub for helpful discussions. Correspondence with Adam Falk and Gino 
Isidori on power corrections are thankfully acknowledged. 
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